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We have introduced a CJT calculation in studying nuclear matter beyond mean field approxi- 
mation. Based on the CJT formalism and using Walecka model, we have derived a set of coupled 
Dyson equations of nucleons and mesons. Neglecting the medium effects of the mesons, the usual 
MFT results could be reproduced. The beyond MFT calculations have been performed by thermo- 
dynamic consistently determining the meson effective masses and solving the coupled gap equations 
for nucleons and mesons. The numerical results for the nucleon and meson effective masses at finite 
^ — ■ temperature and chemical potential in nuclear matter are discussed. 
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I. INTRODUCTION 



Mean field theory (MFT) has an important role in the study of nuclear matter. Based on quantum hadrodynamics 
(QHD) and by mean field approximation, one can obtain the results in good agreement with the properties of the 
nuclear matter, like saturation density and binding energy, etc 0, H, [![. For Walecka model (QHD-I), mean field 
t-H ■ approximation could be regarded as replacing the meson fields by their expectation value Q ■ That is to say treating 
the meson fields as the classical fields while quantum fluctuation ignored. From a Green's function formalism, the 
mean field approximation could also be viewed as resumming all the tad-pole graphs in the nucleon self-energy [y] . If 
one wants to make calculations beyond MFT, one should systematically calculate the higher loop contributions from a 

■ self-consistent resummation scheme. At the same time it is important to preserve the thermodynamic consistency 0]. 
By the one particle irreducible effective action formalism, when evaluated to two loops, the calculations become very 
involved and after taking certain approximations the results are found to be not consistent with the properties of the 
nuclear matter {8]. In principle the higher loop contributions should give reasonable results. Therefore one should find 
systematic and consistent resummation scheme and make proper approximations to obtain soundable results beyond 

■ mean field calculation. 

For the effective action formalism, the effective action r(</>), which is the generating functional of one particle 
irreducible graphs, could been systematically calculated by summing all the relevant Feynman graphs to a given order 
of the loop expansion [t| [l(| HH, 03 • However in this scheme it is not convenient to preserve the consistency of the 
effective action and the full propagators at the same order of the approximation. A generalized effective action for 
composite operators has been developed by Cornwall, Jackiw and Tomboulis (CJT) [l3j]. According to CJT formalism 
the usually effective action T(cj)) is generalized to depend not only on (j>(x), a possible expectation value of the quantum 
field Q(x), but also on G(x,y), a possible expectation value of the time-ordered product T$(x)Q(y). G(x,y) is also 
the full propagator of the field. In this case the effective action T(0, G) is the generating functional of the two particle 
irreducible graphs. An effective potential V(<p, G) can be defined by removing an overall factor of space-time volume 
of the effective action. Physical solutions demand minimization of the effective potential with respect to both (f> and 
G fijl, which means 



and 

dVtf, G) 
dG 



(1) 



0. (2) 



This formalism was originally written at zero temperature. Then it had been extended at finite temperature by 
Amelino-Camelia and Pi where it was used for investigations of the effective potential of the A</> 4 theory [l4| and 
gauge theories 15]. In the study of the spontaneous symmetry breaking system, <fi is also the order parameter of the 
transition and nonzero; while in the system without spontaneous symmetry breaking, this parameter is zero. For the 
thermodynamic system the effective potential could be identified as the thermodynamic potential. One can use the 
CJT formalism to study nuclear matter. In our previous work [l(|, we have already applied the CJT formalism to the 
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study of nuclear matter. At the first step, after the proper approximations we have successfully reproduced the MFT 
results through the CJT formalism. And we have also illustrated that the thermodynamic consistency is preserved in 
the CJT formalism. In this paper we will further study how to make thermodynamic consistent calculations beyond 
MFT in the CJT formalism. 

The organization of the present paper is as follows: In section 2, we formulate the CJT formalism in the study 
of nuclear matter by the Walccka model. The thermodynamic potential and coupled Dyson equations of the full 
propagators are derived at Hartree-Fock approximation. In section 3, we will make the beyond MFT calculations by 
including the medium effects of the mesons in a thermodynamic consistent way and give our numerical results. The 
last section comprises a summary and discussion. 



II. CJT FORMALISM IN NUCLEAR MATTER 



We start from the QHD-I Lagrangian which can be written as 

C = #7^" - m N + g a a - guy^W + \{d^a - m 2 a a 2 ) - -F^F*" + (3) 

where ip, a and uj are the fields of the nucleon, sigma meson and omega meson respectively and F^ — d^uj v — d v uj^. 
According to the Lagrangian we can write down the free propagators of the nucleon, sigma meson and omega meson 
respectively as 

p — TOj\r 

A„(p) = -a- 5 -* , (5) 
p z — mg. 

D^ip) = ^Cr- (6) 
p A — m^j 

We will use the imaginary-time formalism to compute quantities at finite temperature [171 ] . Our notation is 

where (3 is the inverse temperature, (3 = l/T. We have u>„ = 2mtT for boson and w n = (2n + 1)ttT for fermion, where 
n = 0, ±1, ±2, • ■ ■ . A baryon chemical potential \i can be introduced by replacing p Q = iui n with p Q = iu> n + p in the 
nucleon propagator. 

According to the CJT formalism [13] , the expansion of the effective potential or the thermodynamic potential of 
the nuclear matter can be written as 

n(G,A,D) = i f Tr\nG (p)G- 1 (p) + i [ Tr [G^{p)G(p) - l] - % - f In Aq^A^) -if [A^(p)A(p) - l] 
Jp Jp ^ Jp ^ Jp 

- iJ\nD {p)D-\p) -\J[D^{p)D{p) 1] + Q 2 (G,A,D), (8) 

where G, A and D are the full propagators of nucleon, a meson and u> meson respectively, which are determined 
by the stationary condition f2 2 (G, A, D) is given by all the two-particle irreducible vacuum graphs with all the 
propagators treated as the full propagators. In the CJT formalism at Hartree-Fock approximation ^(G, A, D) can 
be illustrated by the graphs of figure [U where the solid lines repent G(p), the dashed line represents A(p) and the 
wavy line represents D tiy {p). The vertices for aip-tp and are — g a and g^y 1 * respectively. The analytic expression 

is 

n 2 (G, A,D) = t -f JjTr [G(q)G(q - p)A(p)] + l A JJtt [y»G{q)y»G{q p)D, v (p)] . (9) 
From the stationary condition @ which demands that f2(G, A, D) be stationary against variations of G, A and D 



q-p q-p 

FIG. 1: Hartree-Fock approximation to Q2(G, A, D). The solid, dash and wavy lines are the nucleon propagator G, the a 
meson propagator A and the u> meson propagator D respectively, q and p are the four momenta. 

respectively, we will have the following Dyson equations 

G(q)- 1 = Goiq^+gl [ G(q~ P )A(p)+gl f ^G{q- p)^ D^ijp), (10) 

J p J p 



A(p)- 1 = Aoip)- 1 -gl Tr [G(q)G(q - p)] , (11) 

Jq 

D^p)- 1 = Do^ip)- 1 - gl f Tr [ %l G(qh„G(q - p)] . (12) 

Jq 

The above equations can be also represented pictorially in figure [2] These integration equations are nonlinear coupled 




D-\p) = D \p) + 



q-p 



FIG. 2: The Dyson equations satisfied by the nucleon, the a meson and the lu meson propagators at Hartree-Fock approximation. 
q and p are the four momenta. 

and momentum dependent. Needless to say they are very difficult for computation. The certain approximations 
should be adopted. 

III. THE BEYOND MEAN FIELD CJT CALCULATION AND NUMERICAL RESULTS 

In our previous study in order to reproduce the mean field results through the CJT formalism, we have 
neglected the medium effects of the mesons which means the meson full propagators are replaced by their bare ones. 
The Dyson equations are decoupled and can be numerically solved. Then we could reproduce the correct saturation 
properties of nuclear matter as shown in figure[3l The coupling constants thus are fixed at g^ = 155.6 and gl — 521.5. 
The result from MFT is also presented for comparison in figure [3] 

In this paper, we are going to make calculations beyond mean field. We will take the following ansatz of the meson 
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FIG. 3: The average energy per nucleon minus the nucleon mass mjv as functions of baryon density at zero temperature in the 
CJT formalism (solid line) and MFT (dashed line). 



full propagators 



D "» = (») 

where M a and M u are the effective masses of the a and w mesons, which in principle should be determined by 
equation (lll| and ()12[) . However, as we have mentioned that the momentum dependent coupled Dyson equations are 
very difficult to be solved. So we have to find other simplified ways to determine the effective masses and meanwhile 
preserve the thermodynamic consistency. 

The full nucleon propagator is taken as the following form 

7 M g" - m N + E 

where E is the proper nucleon self-energy which will be determined by the equation (|10[) . It can be generally written 
as 

E(g) = E.(?) - 7 °Eo(g) + 7-qE„(?)- (16) 
Thus we can define an effective nucleon mass 

M N (q) =m N -E s (g), (17) 

which is momentum dependent. As the nuclear matter is a uniform system at rest, the contribution of E„ term in 
equation (|16p is very small and usually neglected [f|. So only E s and Eo need to be evaluated. 
By the above ansatz of the full propagators, the thermodynamic potential could be written as 

Q(G,A,D) = i JjrlnG (p)G- 1 (p)- Tr [G(p)E(p)] - % - jhx Ao^A" 1 ^) - i jf (Af CT 2 - m 2 ) A(p) 

' / "lnU (p)I>- :l (p) + J [ (M*-ml)g^D^(p)+n 2 (G,A,D). (18) 



2j p — -"2., 
From equation (TTU|) . ([T5]) and (fJ/J]), we could derive that 

E(g)=ifl* I G(q-p)A(p)+ /,-: / -„,<VV/-7>»vI>"». I I.')! 

thus we have 



n 2 (G,A,D) = -TrJ G(q)^g 2 a J G(q - p) A(p) + igi J lfl G(q - p) lv D^ '(p) 

= 5/ Tr[G(g)E( 9 )]. (20) 
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In the next we will mainly calculate E, M„, and in a thermodynamic consistency way. In the following calcu- 
lations the vacuum fluctuations will be neglected. The £(<?) could be evaluated from equation (TT5I) . The momentum 
dependence can be eliminated if we define the nucleon effective mass as the pole mass of the full propagator [l8| in 
the limit q — > as discussed in [lj| . With the £„ term dropped we have 

S = E S - 7 °E , (21) 
_ 2 f d3 P M N 2{M N - E)n a - E a {fi + + n_) 2 f d 3 p 2M N 2(M N - E)n u - E u (n+ + n_) 

{2-Kf2E <J E {M N -Ef-El 9 " J (2nf E^E {M N - Ef - E* A ' 

y = 2 f d3 P -{n+-n-) 2 f d 3 p 

5ff 7 (2tt) 3 2(M N - E f - 2EI +9 "J (2tt) 3 (M n - Ef - E*' [ ' 



where E a / U — Jp 2 + "m 2 , and E = \Jp 2 + M^; the Bose and Fermion distribution functions are 



r/ 

1 _ _ 1 

~ e PE a/ „ _ I ' n ± - e /3(B TAl *) + J' ( 24 ) 

where /i* = — Eo and /x is the baryon chemical potential. 

After performing the frequency sums in equation (|18p and neglecting the vacuum fluctuations, the thermodynamic 
potential could be evaluated as 

n( M) T)= - 4T / -fP\ln(l + e-K E -^)+ln(l + e-K E+ ^) 
J (2ny L 



J "' ~ ' (^^ <s++s -'- 4S »/(^ (fi ' 



4T / W Ml - ."«■) - 2<«* -^SwfK + °* (25) 



where 



The E s and So m ^2 will be evaluated in the forms of equations (|22p and (|23|) . Now we are in a position to determine 
the effective masses M a and M w . In equation (|2l)|) , M CT and M u could be viewed as independent variables. As 
known that a thermodynamic system in equilibrium will minimize its thermodynamic potential to these independent 
variables, which means 

^ 0, 0. (27) 



dM a dM*. 

M a and M u determined in this way are momentum independent and the thermodynamic consistency is certainly 
preserved. From equation (|25j) . the two equations in ([27)) could be written as 



1/71^2 2\ ® f d 3 p n, { nl# /" d 3 p 1 _ ^ <9£ s o f d 3 p ( _ _ ^ 5E 



-(Mf - m^)- / t-^t— + 2M N / — -Vt — (n+ + — - 2 / — -^(n+ - ft-)— - = 0, (28) 

2 V dM a J (2tv) 3 E a J (2ir) 3 E y + ' dM a J (2tt) 3 V + 9M CT ' V ' 



The equations (|2"2")) , (|2"3"|) , (|2"5j) and form a set of closed equations from which E s , So, M CT an d M u could 
be numerically solved. At this stage we could make calculations and discussions beyond MFT in studying nuclear 
matter. As an application we will discuss the temperature and chemical potential dependence of the effective masses 
of nucleons, a mesons and u) mesons in the hot and dense nuclear matter. 

At zero chemical potential, we could discuss the temperature dependence of the effective masses. In figure U the 
nucleon effective mass of the CJT case decreases with temperature increasing. At low temperature the effective mass 
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FIG. 4: The effective nucleon masses as functions of temperature at zero chemical potential in the CJT formalism (solid line) 
and MFT (dashed line). 

decreases very slowly. At certain high temperature the effective mass becomes decreasing fast with temperature 
increasing. The MFT result is also plotted for comparison. We could see that when temperature is high the effective 
mass decreases more rapidly in the CJT case than that of the MFT case. In figure [5] the meson effective masses 
as functions of temperature are plotted. The effective mass of a meson increases with temperature increasing. At 
high temperature it increases more obviously. While The effective mass of lu meson will decrease with temperature 
increasing. It decreases more quickly at high temperature. 
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FIG. 5: The effective masses of a and uj mesons as functions of temperature at zero chemical potential in the CJT formalism. 

At given temperature, we could also study the chemical potential dependence of the effective masses. In figure|6]we 
could see that the nucleon effective masses decrease with chemical potential increasing at different given temperatures. 
From low temperature to high temperature the effective masses decrease more and more remarkably with chemical 
potential increasing. In the CJT case the effective masses decrease more rapidly with chemical potential than that in 
the MFT case at high temperatures. 

We show in figures [7] and [5] the a and to effective masses as functions of chemical potential at different give 
temperatures. The a effective masses increase with the increase of the chemical potential while the to effective masses 
decrease with the increase of the chemical potential at given temperature. For both a and to mesons, the effective 
masses change slightly with chemical potential at low given temperatures but rapidly at high given temperatures. 

IV. SUMMARY 

By the CJT formalism and Walecka model, we have made a beyond MFT calculation in studying nuclear matter. 
We derive a coupled Dyson equations for the full propagators of nucleons and mesons. After neglecting the vacuum 
fluctuations and the momentum dependence in the nucleon self-energies, the nucleon self-energies can be evaluated 
by the Dyson equation of the nucleon full propagator, while the meson effective masses are determined by minimizing 
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FIG. 6: The effective nucleon masses as functions of chemical potential at different given temperatures in the CJT formalism 
(solid line) and MFT (dashed line). 
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FIG. 7: The a effective masses as functions of chemical potential at different given temperatures in the CJT formalism . 



the thermodynamic potential to these effective masses instead of directly evaluating the meson Dyson equations. This 
kind of treatment simplifies the calculation meanwhile preserves the thermodynamic consistency. The effective masses 
and thermodynamic potential can be numerical evaluated at given chemical potential and temperature. The numerical 
results show that the nucleon effective mass decreases with the increase of temperature or chemical potential, and 
it decreases more rapidly in the CJT case than that of the MFT case. For the mesons, the effective mass of oj 
meson decreases with temperature or chemical potential increasing, while the effective mass of a meson increases with 
temperature or chemical potential increasing. Both nucleon and meson effective masses change more quickly at high 
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FIG. 8: The u> effective masses as functions of chemical potential at different given temperatures in the CJT formalism 
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temperature or chemical potential than that at low temperature or chemical potential. 



Acknowledgments 



This work was supported in part by the National Natural Science Foundation of China with No. 10547112 and No. 
10675052. 



[1] R.J. Furnstahl, B.D. Serot, Comments Nucl.Part.Phys. 2 (2000) A23. 
[2] B.D. Serot, J.D. Walecka, Int.J.Mod.Phys. E6 (1997) 515. 
[3] B.D. Serot, Rep.Prog.Phys. 55 (1992) 1855. 

[4] J.D. Walecka, Theoretical Nuclear and Subnuclear Physics, Oxford Univ. Press, 1995. 
[5] B.D. Serot, J.D. Walecka, Adv.Nucl.Phys. 16 (1986) 1. 
[6] S.A. Chin, Ann.Phys. 108 (1977) 301. 
[7] H. Uechi, Phys. Rev. C41 (1990) 744. 

[8] R.J. Furnstahl, R.J. Perry, B.D. Serot, Phys.Rev. C40 (1989) 321. 

[9] R. Jackiw, Phys.Rev. D9 (1974) 1689. 
[10] L. Dolan, R. Jackiw, Phys.Rev. D9 (1974) 3320. 
[11] S. Coleman, E. Weinberg, Phys.Rev. D7 (1973) 1888. 
[12] A. Linde, Rep.Prog.Phys. 42 (1979) 390. 

[13] J. Cornwall, R. Jackiw, E. Tomboulis, Phys.Rev. D10 (1974) 2428. 
[14] G. Amelino-Camelia, So- Young Pi, Phys.Rev. D47 (1993) 2356. 
[15] G. Amelino-Camelia, Phys.Rev. D49 (1994) 2740. 
[16] Song Shu, Jia-Rong Li, Nucl.Phys. A760 (2005) 369. 

[17] J. Kapusta, Finite- Temperature Field Theory, Cambridge Univ. Press, 1989. 
[18] S. Gao, Y.-J. Zhang, R.-K. S, Phys.Rev. C52 (1995) 380. 



